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Abstract 

We show that it is possible to obtain credible static anisotropic spherically symmetric matter 
configurations starting from known density profiles and satisfying a nonlocal equation of 
state. These particular types of equation of state describe, at a given point, the components 
of the corresponding energy-momentum tensor not only as a function at that point, but as a 
functional throughout the enclosed configuration. To establish the physical plausibility of the 
proposed family of solutions satisfying nonlocal equation of state, we study the constraints 
imposed by the junction and energy conditions on these bounded matter distributions. 

We also show that it is possible to obtain physically plausible static anisotropic spher- 
ically symmetric matter configurations, having nonlocal equations of state, concerning the 
particular cases where the radial pressure vanishes and, other where the tangential pres- 
sures vanishes. The later very particular type of relativistic sphere with vanishing tangential 
stresses is inspired by some of the models proposed to describe extremely magnetized neutron 
stars (magnetars) during the transverse quantum collapse. 



1 Introduction 

Classical continuum theories are based on the assumption that the state of a body is de- 
termined entirely by the behavior of an arbitrary infinitesimal neighborhood centered at 
any of its material points. Furthermore, there is also a premise that any small piece of the 
material can serve as a representative of the entire body in its behavior, hence the govern- 
ing balance laws are assumed to be valid for every part of the body, no matter how small. 
Clearly, the influence of the neighborhood on motions of the material points, emerging as 
a result of the interatomic interaction of the rest of the body, is neglected. Moreover, the 
isolation of an arbitrary small part of the body to represent the whole clearly ignores the 
effects of the action of the applied load at a distance. These applied loads are important 
because their transmissions from one part of the body to another, through their common 
boundaries affect the motions, hence the state of the body at every point. The relevance of 
long-range or nonlocal outcomes on the mechanical properties of materials are well known. 
The main ideas of Nonlocal Continuum Mechanics were introduced during the 1960s (see |j] 
and references therein) and are based on considering the stress to be a function of the mean 
of the strain from a certain representative volume of the material centered at that point. 
Since then, there have been many situations of common occurrence wherein nonlocal effects 
seem to dominate the macroscopic behavior of matter. Interesting problems coming from a 
wide variety of areas such as damage and cracking analysis of materials, surface phenomena 
between two liquids or two phases, mechanics of liquid crystals, blood flow, dynamics of col- 
loidal suspensions seem to demand this type of nonlocal approach which has made this area 
very active concerning recent developments in material and fluid science and engineering. 

The true equation of state that describes the properties of matter at densities higher 
than nuclear (« 10 14 gr/cm 3 ) is essentially unknown due to our inability to verify the 
microphysics of nuclear matter at such high densities @-|[J. Currently, what is known 
in this active field comes from the experimental insight and extrapolations from the ultra 
high energy accelerators and cosmic physics (see j?| and references therein). Having this 
uncertainty in mind, it seems reasonable to explore what is allowed by the laws of physics, 
in particular considering spherical or axial symmetries within the framework of the theory 
of General Relativity. In the present paper, we shall consider a Nonlocal Equation of State 
(NLES from now on) in general relativistic static spheres. It follows a previous work || on 
collapsing configurations having a NLES. The static limit of the particular NLES considered 
in this reference can be written as 
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p r (r) = p (r) - -^ / r 2 p(r) dr + -^— ; (1) 



where C is an arbitrary integration constant. It is clear that in equation ([!]) a collective 
behavior on the physical variables p(r) and P r (r) is present. The radial pressure P r (r) is not 
only a function of the energy density, p(r), at that point but also its functional throughout 
the rest of the configuration. Any change in the radial pressure takes into account the effects 
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of the variations of the energy density within the entire volume. 

Additional physical insight of the nonlocality for this particular equation of state can be 
gained by considering equation (D) re-written as 



*W = *)-§<,>, + £. with Wr = £±^=*M. (2) 



Clearly the nonlocal term represents an average of the function p(r) over the volume enclosed 
by the radius r. Moreover, equation (§) can be easily rearranged as 



12 C 1 2 C 

gP(r) + 3 (p(r) - (p(r)» + ^ = gP(r) + - o p + ^ 



^r(r) = -p(r) + - {p{r) - (p(r)» + -3 = -p(r) + - ^ + -3 , (3) 



where we have used the concept of statistical standard deviation o p from the local value of 
energy density Furthermore, we may write: 

( J>(r) = lp(r) and 
P r (r) =V(r)+ 2a v{r) + — where I 

{ a T{r) = (|p(r) - I (p) r ) = (V(r) - P(r)) . 

(4) 

Therefore, if at a particular point within the distribution the value of the density, p(r), gets 
very close to its average (p{r)) the equation of state of the material becomes similar to the 
typical radiation dominated environment, P r (r) ~ V(r) = \p{r). 

In reference 0, it is shown that under particular circumstances a general relativistic 
spherically symmetric anisotropic (nonequal radial and tangential pressures, i.e. P r 7^ P±) 
distribution of matter could satisfy a NLES. Some of these dynamic bounded matter configu- 
rations having a NLES with constant gravitational potentials at the surface, admit a Confor- 
mal Killing Vector and fulfill the energy conditions for anisotropic imperfect fluids. Several 
analytical and numerical models for collapsing radiating anisotropic spheres in general rel- 
ativity were also developed in that paper. Although the perfect pascalian fluid assumption 
(i.e. P r = P±) is supported by solid observational and theoretical grounds, an increasing 
amount of theoretical evidence strongly suggests that a variety of very interesting physical 
phenomena may take place giving rise to local anisotropy. In the Newtonian regime, the 
consequences of local anisotropy originated by anisotropic velocity distributions have been 
pointed out in the classical paper by J.H. Jeans 0. In the context of General Relativity, 
it was early remarked by G. Lemaitre [10 1 that local anisotropy can relax the upper limits 



imposed on the maximum value of the surface gravitational potential. Since the pioneering 
work of R. Bowers and E. Liang |l]J , the influence of local anisotropy in General Relativity 
has been extensively studied (see |T2 and references therein). Recently, the role of local 



anisotropy in the equation of state modeling the interior of extremely magnetized neutron 
stars or Magnetars at supranuclear densities, has been attracting the attention of researchers 



(see ||13||-fl7] and references therein). It has been shown that these intense magnetic fields 
could induce local anisotropy in pressures within the equation of state of the matter dis- 
tribution. In fact, two of these references ||14|| and ]16[ consider the consequences on the 



instability of the matter configuration in the most extreme case of anisotropy in pressures, 
i.e. when there exists a vanishing pressure perpendicular to the magnetic field. In this case 
the vanishing of the pressure transverse to the magnetic field leads the configuration to a 
transverse quantum collapse of the configuration. 

The present paper is focused on the plausibility of obtaining static solutions satisfying 
NLES and whether those solutions could represent a credible bounded anisotropic matter dis- 
tributions in General Relativity. We shall show that it is possible to obtain static anisotropic 
spherically symmetric matter configurations starting from known density profiles and satis- 
fying a NLES. We also present two static solutions for NLES configurations concerning the 
particular case of vanishing radial pressure (P± 7^ and P r = 0) and vanishing tangential 
stresses, i.e., Pj_ = and P r ^ 0. 

The structure of the present work is the following. The next section contains an outline 
of the general conventions, notation used, the metric and the corresponding field equations. 
Section |2| is devoted to solve the Einstein Field Equations for a matter distribution assuming 
a nonlocal equation of state. The family of NLES is presented in Section ^. In Section ^ we 
study the consequences imposed by the junction and energy conditions on bounded matter 
distribution. The method and some fluid sphere models are considered in the Section [5[ The 
two particular cases, i.e. P r = 0, and P± = 0, are presented in Section || Finally, in the last 
section, our concluding remarks and results are summarized. 

2 The Einstein Field Equations 

To explore the feasibility of nonlocal equations of state for bounded configurations in General 
Relativity, we shall consider a static spherically symmetric anisotropic distribution of matter 
with an energy-momentum represented by T^ = diag (p, — P r , —P±, —P±), where, p is the 
energy density, P r the radial pressure and P± the tangential pressure. 

We adopt standard Schwarzschild coordinates (t, r, 8, <p) where the line element can be 
written as 

ds 2 = e 2 ^dt 2 - e 2X ^dr 2 - r 2 dfi 2 , (5) 

with dfl 2 = d9 2 + sin 2 9d(f) 2 , the solid angle. 



The resulting Einstein equations are: 
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where primes denote differentiation with respect to r. 

Using equations (0) and @, or equivalently the conservation law T£. 
the hydrostatic equilibrium equation for anisotropic fluids 



0, we obtain 



P' r = -(p + P r )u' + - (P ± - P r 
r 

Equation (Rl) can be formally integrated to give 



(9) 
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m(r) 
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where a mass function m(r), has been defined by 



m(r 



An I pr 2 dr 
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and it corresponds to the mass inside a sphere of radius r as seen by a distant observer. 

Finally, from (|9|), ( |10D and (^) the anisotropic Tolman-Oppenheimer-Volkov (TOV) equa- 
tion IjT^I can be written as 



dP r , _. /m + 4vrr 3 P r \ 2 /r> 

dr \r [r — 2m) J r 



(12) 



Obviously, in the isotropic case (Pj_ = P r ) it becomes the usual TOV equation. 

It has been established |[8||-j2(| that if p is a continuous positive function, Pj_(r) is 
a continuous differentiate function and P r {r) is a solution to the equation with starting 
value Pj_(0) = P r (0), there exists a unique global solution to (O) representing a spherically 
symmetric fluid ball in General Relativity. 



3 A Family of Solutions with a NLES 

In this section we are going to present a family of static solution of the Einstein equations 
satisfying a NLES. Defining the new variables: 



^=Mr)e^, and e 2A « 



h (r) 



with h(r) = \-2 



m (r) 



the above metric (|5T) can be re-written as 



(13) 



ds 2 = h (r) e A ^dt 2 - ^-dr 2 - r W 



h(r) 



The resulting Einstein Equations are: 
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Now, equation ([J) can re-stated as a differential equation 
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r"p{r) dr + - 
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P r (r) = p(r) 3 / f 2 p{f) df + 

r Jo 
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<* p-3P r + r(p' - P' r ) = 



Thus, from (|18l), (|X5|) and (0), we have 



- (/i' + 2h{3') + //' + 2/3W + 2h{3" = 
r 



which can be formally integrated yielding 

l3(r) 



M£' 
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dr + Ci 



where C and C\ are arbitrary integration constants. 

At this point, equation (^) deserves several comments: 

• Firstly, if we set C = C\ = , a particular family of solutions 



(17) 



(19) 



(20) 
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considered in a previous work || is found. The metric obtained in this particular 
static case recalls the so called isothermal coordinates system [EH]] which in turn is a 
particular case of the more general "warped space-time" (we refer the reader to p2 
and references cited therein for a general discussion) 

Next, is the approach we have followed in order to obtain static anisotropic solutions 
having a NLES. It is clear that if the profile of the energy density, p(r), is provided, the 
metric elements h(r) and j3(r) can be calculated through (|TT|), dl3|) and (|20|) . Therefore, 
we can develop a consistent method to obtain static solutions having NLES from a 
known static ones. 



Finally, the metric elements flT3|) and fl20|) , should fulfill the junction conditions and 
the physical variables coming from the energy momentum tensor, are only restricted by 
the hydrostatic equilibrium equation flTJj) and by some elementary criteria of physical 
acceptability. The next section is devoted to list these criteria for anisotropic and 
isotropic fluids with a NLES. 



In terms of the metric elements ([13]) and |2(J, the corresponding Einstein Field Equations 
for anisotropic fluids having NLES can be written as: 



inp 



2m! 



(22) 



in 



Pr 



2m! 4 (m - C) 



and 



(23) 



„ m" 2 (m'r — m) 

8nP ± = — + -i >- 



m r 



m 



r — 2m 



(24) 



4 Junction and Energy Conditions 



Most exact solutions of the differential equation ( |12D supplied by the literature have been 
obtained from excessively simplifying assumptions solely with the purpose to find such solu- 
tions, and, consequently, they rarely represent physically "realistic" fluids (see for example, 
two interesting and complementary reviews on this subject: |23| and [^4|). In order to es- 
tablish the physical acceptability of the proposed static family of solutions fll~3|) and (|20|). 
satisfying a NLES, ([!]) we state the consequences imposed by the junction and energy con- 
ditions for anisotropic fluids on bounded matter distribution. 

We are going to consider bounded configurations, i.e. a matter distribution isolated in 
the sense that, where the pressure vanishes occurs at a finite radius. In fact, it can easily 



be shown that the necessary and sufficient condition for matching an interior solution (|i~4"D 
onto the exterior Schwarzschild solution of total mass M, 



ds 2 
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M 



dt 2 - 1 



M\ _1 , 9 
2— dr 2 
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(25) 



is that the pressure equals zero at a finite radius r = a. In this case, the interior metric ([14]) 
should satisfy the following conditions at the boundary surface of the sphere r = a : 

(26) 



P r (a) = => p(a)=/3 a = 0, and 
^From the condition P r (a) = 0, it is clear that 

M = 4vr / p f 2 df = 2vra 3 p (a) . 
Jo 

Now, equation (|23|) leads to C = i.e. 

m(r) =C + 2vr[p(r) -P r (r)]r 3 if m(0) = 
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We should require that the metric elements have to be finite and non zero everywhere 
within the matter configuration, with no changes of sign nor loss of reality allowed. This 
implies that 



m(r) > - 
v ; 2 



V r 



(29) 



where m(r) is the mass function defined by equation (|i~l~D . Recently |2(J the above condition 
(|9|), for spherically static bounded configurations, has been proved to be equivalent to 
requiring that p + P r + 2 P± > 0, (i.e. the strong energy condition ) which clearly 
implies that the space time cannot contain a black hole region. 

Next, we will establish the consequences of the physical acceptability for an anisotropic 
fluid configuration in terms of the mass function m(r). 

1. The density must be positive definite and its gradient must be negative everywhere 
within the matter distribution. Trivially, from equation (p2[) and its derivative we 
obtain 

on 2ttz 
p > => rri > and -£- < => m" < . (30) 



2. The radial and tangential pressure must be positive definite. Therefore, equations (|23| 
and pj), yield 
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3. The radial pressure gradient must be negative definite. Thus, differentiating equation 
we get 

dP r „ „ 4m' 6m 

IT " 



< 



m" < 



(33) 



4. The speed of sound should be subluminal, consequently from equation (|18|) we obtain 

- dPr .. (P"3P r ) _ 



f.. = 



dp'' * r(£) 
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It is clear that, due to the density gradient being negative everywhere within the 
configuration, the requirement of subluminal sound speeds leads to p > 3 P r for both 
anisotropic and isotropic fluids having a NLES. 

In addition to the above intuitive conditions we should satisfy either the Strong Energy 
Condition or Dominant Energy Condition: 

• Demanding that the trace of the energy-momentum tensor be positive we find for the 



Strong Energy Condition, p + P r + 2P_ L > 0. Thus, using equations (p2|), (p3|) and (^4| 
we get 



p + P r + 2 P ± > 
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The Strong Energy Condition also implies p + P r > and p + P± > 0, therefore, 
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In the case of isotropic fluids, the above condition ( J35|) reads 
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(38) 



The Dominant Energy Condition entails that the density must be larger than the 
pressure. Now, by using equations fl2"2|), (|2"3"| ) and (^) we get 
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Notice that, if the density, the radial and the tangential pressures are positive definite 
and the Strong Energy Condition is satisfied, then the Dominant Energy Conditions ( P§D 
and fllOD are automatically fulfilled, but the inverse is not true. 

In order to determine the physical reasonableness of the anisotropic configurations having 
a NLES we shall explore two different sets of conditions. Both sets include conditions ( p9[) 
through (J34[) but differ in the selection of the Strong Energy Condition or the Dominant 
Energy Condition. In the first group, we shall use the Strong Energy Condition fl35|), and 
because we require subluminal sound velocities (@), p > 3 P r will be required. Obviously, 
when isotropic fluids are considered, the Strong Energy Condition reads p > 3P r . and the 
sound speed will be subluminal in all cases. The second set again comprises conditions ( f29|) 
through fl3"4|), the fulfillment of the Dominant Energy Conditions (equations fl39| ) and (f40|)) 
and the requirement of p > 3P r . 

5 A Method for NLES Static Solutions 

In the present section we shall state a general method to obtain NLES static anisotropic 
spherically symmetric solutions from known density profiles. 

5.1 A Method for NLES Anisotropic Solutions 

Concerning anisotropic fluids, i.e. non pascalian fluids where P r ^ P±, the method to obtain 
NLES static spherically symmetric solutions is as follows: 

1. Select a static density profile p(r), from a known static solution. Then, the mass 
distribution function, m = m(r), can be obtained through equation (pTTJ) . The junction 
condition (^) implies the continuity of h(a) and the expression for the total mass 
m(a) = M can be procured. 

2. Next, check out where and under what circumstances, all the above physical and energy 
conditions, written in terms of the mass function, m = m(r), are fulfilled. In other 
words, the mass distribution function obtained from the density profile selected should 
satisfy the inequalities fl2jj) thought (|5B|) (or (f2|) thought (g§ and (§g) and fl^)), at 
least for some region [ri,^] within the matter distribution with < r\ < T2 < a, and 
for particular values of the physical parameters that characterize the configuration. 



3. Following, the other metric coefficient, /3(r), can be found by using equation (|20|). 
Notice that because the boundary conditions (|28|) the actual expression for f3(r) will 
be (©. 



4. Finally, Einstein Field Equations (|23| ) through fl24|) provide the expressions for the 
radial and tangential pressures, P r and P±, respectively. 



5. The integration constants are obtained as consequences of the junction conditions at 
the boundary surface, r = a, i.e. /3(a) = 0, m(a) = M and P r (a) = 0. 

In order to illustrate the above procedure we shall work out several examples for six static 
density profiles borrowed from the literature. 



5.2 Examples of NLES Static Solutions 

Example 1: The first example comes from a density profile proposed by B. W. Stewart 
J25| , to describe anisotropic conformally flat static bounded configurations: 
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with K = const. 

Thus, radial and tangential pressures can be written as 
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The constant K can be found from the boundary condition M = m(r = a) as 



K = — In 

2a 



1 + m^ 

l-(f) 1 



while because the pressure vanishes at the surface P r (r = a) = we have: 
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8Kae 



2Ka 
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(41) 



(42) 



(43) 



(44) 



Example 2: The density profile of the second example is found due to P.S. Florides [E6 



but also, corresponding to different solutions, by Stewart |£5| and more recently by M. 
K. Gokhroo and A. L. Mehra ||27|| . The Gokhroo-Mehra solution, represents densities and 
pressures which, under particular circumstances [pgfl , give rise to an equation of state similar 
to the Bethe-Borner-Sato newtonian equation of state for nuclear matter P, |3|, |29 . 
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(45) 



with a and K = const. 

The radial and tangential pressures are obtained from Einstein Field Equations 
through ( p4|) and can be expressed as 
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(47) 



where a = 0.631515. 

The boundary conditions at the surface r = a, i.e. M = m(a) and P r (a) = 0, lead to 



5 , , , era 3 

K = - and M = 

9 9 



(48) 



Example 3: This solution was discovered by H.B. Buchdahl |$0] and rediscovered later by 
IJ. The corresponding density profile can be written as. 
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(49) 



with C = const. 

The method applied to the above density profile leads to 
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and P\ 
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16tt (2 - Cr 2 ) (1 + Cr 2 ) 



3 • 



(50) 



where C = -^ coming from P r (a) = 



Example 4 The matter distribution sketched in this example is borrowed from Tolman IV 



isotropic static solution which was originally presented by R.C Tolman in 1939 ||32| . Tolman 
IV static solution is, in some aspects, similar to the equation of state for a Fermi gas in cases 
of intermediate central densities. This same profile (and Tolman IV solution) is also found 



as a particular case of a more general family of solution in P3J and 
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Einstein Field Equations (^) 
pressures as 



P 



provide the expressions for the radial and tangential 
C l-2x-K(l + x + 2x 2 ) 
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(1 + 2xf (1 + x) (1 + Kx) 



(53) 
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where x = Cr 2 ;the constants K and C are obtained from the boundary conditions M 
and P r (a) = respectively, i.e. 



K 



xJ-2MC2(l + 2x 1 ] 

3 

x\ (1 + Xi) 



and M 
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Chl + xi + 2x\ 



2 ' 



m(a) 



(54) 



and x\ = C a 2 



Example 5: The density profile of this example corresponds to a solution originally pro- 
posed M. Wyman f3~5fl . Again, the same solution is found in |33|, |34], |36], [37| and |38| 
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1 Kx^ 
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(55) 



with x = C r 2 and K,C = const. 

The method applied to the density profile (|55| ) leads to radial and tangential pressures 
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8 ^(1 + 3^)3 



and 



(56) 



where 
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(58) 



The constant K has been obtained for boundary conditions and the pressure at the surface 
(P r {a) = 0) determines the next integration constant: C := \ 



Example 6: The following density profile was found by M.P. Korkina, in 1981, [Q and 
rediscovered a year later by M.C. Durgapal fl3"4 |. It can be written as 
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Again, for the density profile ([59]) the radial and tangential pressures are found to be 

C 3 (1 + Ax)^ (1 - x) - 2K (1 - x - 8x 2 ) 



8tt 



1 + xf (1 + Ax) 



and 



(60) 



Pi 



C AK 2 xA + AK(l + AxY> 2 + 3x11 + 6 



167T 



l + x) 3 (l+4x)2 (1 + Ax)* (x - 2) - 2Kx 



(61) 



where 



A = -8x 4 + 3Ax 3 + 36x 2 + 13x + 1 



5X 



65x 



16x + 5x + x — 1 



and 



IT = 192x 4 + 80a; 3 + 116x 2 + 87x + 23 . 
Finally, the constants K and C can be obtained for boundary conditions: 

3 x\ 



K 




(l + 4xi)3 




2x\ 


with 


X\ 


= Ca 2 



AMC*(l + xi)-3xf 



and 



M 



C* (8x 2 + x x - I) 



(62) 



5.3 Modeling Anisotropic Spheres with NLES 

The parameters: mass, M,in terms of solar mass M Q , M/a gravitational potential at the 
surface, boundary redshift z a , surface density p a and central density p c , that characterize 
these bounded configurations are summarized in the following table 



Equation of State 


M/a 


M(M Q ) 


Z a 


p a x 10 14 (gr.cm 3 ) 


p c xlO 15 (gr.cm 3 ) 


Example 1 


0.32 


2.15 


0.6 


6.80 


1.91 


Example 2 


0.40 


2.80 


1.2 


8.84 


1.99 


Example 3 


0.38 


2.54 


1.0 


8.04 


2.41 


Example 4 


0.25 


1.69 


0.4 


5.36 


2.00 


Example 5 


0.38 


2.54 


1.0 


8.04 


3.04 


Example 6 


0.35 


2.37 


0.8 


7.75 


2.11 



All these parameters have been tuned up in order to that the mass function satisfies the 



physical and energy conditions, i.e., the inequalities ( p9|) thought (|36"D (or (|29|) thought fl3~4] 
and (|39|) and fl40"|)), for a typical compact objects of radius of a = 10 Km.. 
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6 P r = 0, and P± = static solutions: 

In this section we shall show that it is possible to obtain credible anisotropic spherically 
matter configuration with a NLES for two particular cases: one where the radial pressure 
vanishes, i.e., P± 7^ and P r = and other where the tangential pressures vanishes, i.e., 
P± = and P r ^ 0. 



T0fl , A. Einstein 1^] and P.S. Florides[26 



6.1 P r = and P ± ^ 

The study of matter configurations with vanishing radial stresses traces back to G. Lemaitre 

Non static models have been 

mmm 



in the static case, 
considered in the past in references [TT2| , |4"0| , [4T| , [42"|] and recently in |14|, |40|, jig, |48| concerning 
their relation with naked singularities. More recently conformally flat models with vanishing 
radial pressures has been considered in f43| . We shall present here one of these static solution 
satisfying an NLES. 

From equation (^) it is trivial to see that 

2m' Am 
Thus, density and tangential pressure profiles can be calculated as 



m (r) = Cr 2 



(63) 



P 



C 



and 



P, 



C 2 



2vrr " ' "" 4vr(l - 2Cr) ' 

The constant C can be obtained from the boundary conditions, i.e. 



(64) 



m[a) 



M 



*=4 



(65) 



6.2 P ± = and P r ^ 



The second static solution with a NLES we shall consider is the case Pj_ = and P r ^ 0. 
The rationale behind this supposition is the considerable effort that has been directed in 
recent years to study the effects of intense magnetic fields (B > 10 15 G) on highly compact 
astrophysical objects (|[L3||-|r7| and references therein). In fact, some observations seems to 



confirm that newly born neutron stars with very large surface magnetic fields (magnetar) 
could represent either soft gamma repeaters or anomalous X Ray pulsars [[0J, or both. 
One of the most striking effects of these large magnetic fields on the equation of state for 
superdense matter distributions is that they induce local anisotropic pressures |14|, |16| , jITf and 
it is possible to obtain extreme cases were the pressure perpendicular to the magnetic field 
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could vanish. In our case we shall suppose that tangential stresses vanish in all directions. 
In this sense we consider our solution inspired by the magnetar models 
It is clear that if P± = equation (|4|) leads to 



m" 2 (m'r — 

— + - 3" 


m) 


m'r — m 


-1 


r — 2m 







(66) 



which can be integrated yielding 



r 



I r 

m = - [1 — e 



-2(dr+C 2 )l _ (g^ 

The density and the pressure profiles emerge from equations (0) and ([!]), respectively as: 



p -2(C\r+C 2 ) 

p = — [2rd - 1 + e 2 ^ r+c ^] and 

8nr 2 L J 

„-2(Cir+C 2 ) 

P, = - 2 [2rd + 1 - e 2 ^ r+c ^} . 
8irr z L J 



The constants C\ and C<i can be obtained for boundary conditions at r = a as 

M 



m(a) = M 
PJa) = 



d 



Co 



a (a- 2M) 

-M 1 

+ -In 



and 



a-2M 2 \a-2M 
It is worth mentioning that from equation (^), we have 



(68) 
(69) 

(70) 
(71) 



r = 



m = 



i[2+te-2)ln(l-2f)] 



(72) 



Therefore, matter configurations with NLES having vanishing tangential pressures are only 
possible within a region of the sphere, i.e. 



a 
2 



2 + ( s -2jl^l-2^) 



< r < a, 



(73) 



because, within the range < r < | [2 + (-^- — 2) In (l — 2^)] the mass function becomes 
negative. 
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7 Concluding Remarks 

We presented a method to obtain NLES static anisotropic spherically symmetric exact solu- 
tions starting from known density profiles. It is clear that, when such a density profile, p(r), 
is provided the radial pressure, P r (r), can be obtained from the NLES ([I]) and the tangential 
pressure P±(r) can be solved algebraically from the anisotropic Tolman-Oppenheimer-Volkov 
(TOV) differential equation (|12|) . Althought, we supposed to have found a more general fam- 
ily of solutions for the Einstein Equations, fl20]), than the one considered in || but, due to 
the boundary conditions (|28|) , the actual family of solutions is indeed the previous ( |2lD 
considered in reference (pl|). 

It is also obtained credible anisotropic spherically matter configurations with a NLES for 
two particular cases: one where the radial pressure vanishes, i.e., P± ^ and P r = and 
other where the tangential pressures vanishes, i.e., P± = and P r ^ 0. 

We worked out in detail the junction and energy conditions for anisotropic (and isotropic) 



bounded matter distributions, i.e. fl29|) thought ([40]), establishing their consequences in 
terms of the mass function m(r) for bounded configurations having NLES. In the early days 
of Relativistic Astrophysics, some of these results led to the first general theorems due to 



H. A. Buchdahl[[3(J, |49| and H. Bondi|)(], |51| concerning inequalities limiting the behavior of 
the mass function for compact objects. 

Because the mass function m(r), presented in equations ([4TD through fl59|) and the corre- 
sponding pressure profiles are regular at the origin, the solutions to the Einstein Equations 
are considered unique in the sense stated by Rendall and collaborators [18], [l9| . That is to 
say: for a given value of the central pressure, P = Pj_(r = 0) = P r {r = 0), there exists a 
unique global solution P r (r) to the anisotropic TOV equation (Tl2|) representing a spherically 
symmetric fluid sphere in General Relativity. 

It is clear that the last term in the anisotropic TOV equation (|T2|), ^ (P± — P r ) = A, 
represents a "force" due to the local anisotropy. This "force" is directed outward when 
Pj_ > P r 43> A > and inward if P± < P r -&■ A < 0. As it is apparent from the figure 
and also from the above table, those models with a (average) repulsive "anisotropic" force 
(A > 0) allows the construction of more massive distributions. This picture is more evident 
for solutions where P± ^ and P r = and P± = and P r ^ 0, which are the extreme 
situations. Considering plausible static matter configurations (it is to say those that satisfy 
conditions (^) thought (f40|)) we have 



• more massive configurations if the radial pressure vanishes i.e. A > has a maximum, 

• less massive configurations if the tangential pressures vanishes i.e. A < has a mini- 
mum. 

Notice that Model 4 (NLES Tolman IV anisotropic static solution) and Model 5 (NLES 
anisotropic Wyman like static solution) have "soft" (P± < P r ) cores and "hard" (P± > P r ) 
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outer mantles. It is worth mentioning that models presented here are very sensitive to the 
change of values of the parameters sketched in the table. 

As it can be appreciated from equation (J72[), fluids satisfying NLES and having vanish- 
ing tangential pressures are only possible within a region of the matter configuration with 
| [2 + (-g — 2) In (l — 2^)1 < r < a. Within the inner core, the mass function becomes 
negative. This situation emerges from a limitation of the definition of the Schwarzschild 
mass function (|TT|) when it is apply to anisotropic fluid configurations. It is interesting to 
say that the expression of the Tolman Whittaker mass|52] 

pr 

m TW ( r ) = 4vr / r 2 e (l/+A)/2 (p-P r - 2P ± ) dr = e^ +x ^ 2 (m + 4vrP r r 3 ) (74) 

Jo 

for this static NLES model reads 

m TW (r) = d r 2 (75) 

and precludes non physical situations at the core of the distribution. It is clear that for this 
case the Tolman- Whittaker mass function, itltwi is a more suitable concept of mass than 
the Schwarzschild mass function, m (see []12[ for a comprehensive discussion on this point, 



plentiful of examples). More over, from equations (^) and fl75|) , the expression of the rriTw 
for this particular static NLES model with vanishing tangential stresses is the same that m 
for the model with vanishing radial pressure. 

It can also be mentioned as a curiosity that the relation among the areas under the 
curves displayed in the figure, i.e. f Q — (P± — P r ) , corresponds to the relation among de 
total masses. The more massive the model is the greater the area under the curve - (P± — P r ) 
is. 

In principle, a NLES in isotropic pascalian fluid (P± = P r ) can not be ruled out, but we 
have integrated numerically the standard TOV equation (|T2|), which can be written in terms 
of the mass distribution function as 



, 2m' Am 2 (m'r — m) 
m = 



17-1 iriZ qr*£ 



m'r — m 



r — 2m 




with <^ (76) 



and the inequalities (|29|) thought ( J36|) (or ( p9|) thought ( |34"D and ( |39D & fllPp ) corresponding 
the physical and energy conditions were not fulfilled for any configuration having a NLES. 
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Figure caption 

Figure 1: - (P± — P r ) = A, as function of r/a which represents a "force" due to the local 
anisotropy. This "force" is directed outward when P± > P r <^> A > and inward if 
P± < P r <^> A < 0. Model 4 (NLES Tolman IV anisotropic static solution) and Model 
5 (NLES anisotropic Wyman like static solution) have "soft" (P± < P r ) cores and "hard" 
(P± > P r ) outer mantles. 
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